Pairing in a two-component ultracold Fermi gas: phases with broken space symmetries 
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We explore the phase diagram of a two-component ultracold atomic Fermi gas interacting with 
zero-range forces in the limit of weak coupling. We focus on the dependence of the pairing gap and 
the free energy on the variations in the number densities of the two species while the total density 
of the system is held fixed. As the density asymmetry is increased, the system exhibits a transition 
from a homogenous Bardeen-Cooper-Schrieffer (BCS) phase to phases with spontaneously broken 
global space symmetries. One such realization is the deformed Fermi surface superfluidity (DFS) 
which exploits the possibility of deforming the Fermi surfaces of the species into ellipsoidal form at 
zero total momentum of Cooper pairs. The critical asymmetries at which the transition from DFS to 
the unpaired state occurs are larger than those for the BCS phase. In this precritical region the DFS 
phase lowers the pairing energy of the asymmetric BCS state. We compare quantitatively the DFS 
phase to another realization of superconducting phases with broken translational symmetry: the 
single-plane-wave Larkin-Ovchinnikov-Fulde-Ferrell phase, which is characterized by a nonvanishing 
center-of-mass momentum of the Cooper pairs. The possibility of the detection of the DFS phase in 
the time-of-flight experiments is discussed and quantified for the case of 6 Li atoms trapped in two 
different hyperfine states. 



I. INTRODUCTION 

The progress achieved in recent years in trapping and 
manipulating ultracold fermion gases has focused much 
theoretical attention on the pairing properties of di- 
lute fermionic systems. Current experiments with 6 Li 
and K atoms are carried out at temperatures which 
are a fraction (~ 0.1 — 0.3) of the Fermi tempera- 
ture 

[IBSSlHllSli These systems are thus char- 
acterized by a filled Fermi sea, and at sufficiently low 
temperatures attractive two-body interactions are ex- 
pected to drive the Cooper pairing instability 0,0. The 
strength of the two-body interactions can be tuned using 
the Feshbach resonance mechanism by varying the exter- 
nal magnetic field |TlllT3,IT^ : thus, the entire range from 
weak to strong couplings can be probed. In the crossover 
region the Feshbach resonance may strongly enhance the 
pairing interaction and give rise to high-temperature su- 
perfluidity El El El 111 El El MM El El El El 

Recent experiments have probed the condensation of 
fermionic pairs above the Feshbach resonance, where 
the system does not support a genuine two-body bound 
state El El El- The nature of these pairs (tightly 
bound molecules versus extended Cooper pairs) is not 
clear yet, and the signatures of the superfluid phase tran- 
sition, which are manifest in the collective and hydro- 
dynamic behavior of these systems j^, E2' have been 
searched for. The measured collective modes of 6 Li were 
found to be consistent with superfluid hydrodynamics 
and provide some evidence for superfluidity in a reso- 
nantly interacting Fermi gas |28j ] . A gap in the quasipar- 
ticle spectrum of a two-component gas of 6 Li was recently 
observed using radio- frequency spectroscopy [33|. 

The s-wave interaction dominates the pairing inter- 
action in cases where (i) the pairing is among atoms 
of equal mass number but belonging to different hyper- 



fine states or (ii) the pairing is between different species 
(in mixtures of atoms- e. g,. 6 Li and 40 K), in which 
case the symmetric BCS limit is not realized 34]. Sys- 
tems where two hyperfine levels are populated have been 
created and studied experiment ally with 6 Li and 40 K 

atoms EII1I1I1I1I10,I1E1EEE1E1E3 These 
systems are characterized by a hierarchy of scales: the 
typical range of the van der Waals forces is R < 10~ 6 
cm while the de Broglie wavelength of particles at the 
top of the Fermi sea is fcp ~ 10 3 -10 4 cm' 1 . Since 
kpR <C 1 the interaction can be approximated by a zero- 
range force which is characterized by the s-wave scatter- 
ing length as- We specify our discussion below to the case 
where two hyperfine states of 6 Li are populated, in which 
case the scattering length in units of the Bohr radius is 
as/as = —2160. For typical values of Fermi momenta 
quoted above kpas — 0.04 and the system is in the weak- 
coupling regime, since v(kp)Uo — (2/7r)fcF|as| <C 1, 
where v(kp) = mkp / '(2ir 2 h 2 ) is the density of states at 
the Fermi surface and Uq — Airtfas/m is the strength of 
the contact interaction. For larger values of kpas ~ 1 
the weight of the negative-energy states to the single- 
particle spectral function is not negligible - the bound 
states need to be incorporated in the theory along with 
the pair correlations on the same footing [3q|. 

The BCS theory predicts a suppression of the pair- 
ing correlations when the Fermi energies or, equivalently, 
the densities of the two hyperfine states |1) and |2), de- 
noted below as p-\n, are different. In the low-density and 
weak-coupling limit (fci?|as| -C 1) the value of the critical 
density asymmetry a = (p-j — Px)/(PT +Pi)> f° r which the 
superfluidity vanishes, follows from the relation [HE3 
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where A — 8e 2 pexp [— n / (2kp\as\)] is the gap in the 
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symmetric system and p is the chemical potential. There- 
fore, the gap disappears for asymmetries a > a^ cs — 
3Ao/(4/i). For example, if the pairing is between 6 Li 
atoms in the states |1) = \F = 3/2, rap = 3/2), and 
|2) = 1 3/2 , 1/2), in which case the pairing interaction is 
characterized by a scattering length as/ as — —2160, the 
maximum asymmetry at which BCS pairing is possible 
is af cs ~ 0.055 for a density p = p T + p l = 3.8 x 10 12 
cm -3 (corresponding to kp\as\ — 0.55). 

One purpose of this work is to show that the pair- 
ing correlations in ultracold atomic gases can persist for 
density asymmetries a > cef cs and can be enhanced for 
a < a^ cs if the Fermi spheres of two hyperfine states 
are deformed into ellipsoids in momentum space. In the 
strong-coupling regime (Ao//i ~ 0.3) the superconduct- 
ing state featuring deformed Fermi surfaces was found 
preferable to the spherically symmetric BCS state |38| . 
Here we rather explore the weak-coupling regime. Alter- 
natively, the homogeneous BCS phase can evolve into the 
Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) phase [MEI, 
which again sustains a > a^ cs asymmetries by allowing 
the Cooper pairs to carry finite center-of-mass momen- 
tum (the Fermi surfaces in the LOFF phase are spheri- 
cal). 

There has been much interest in the LOFF phase 
outside the condensed matter context in connection 
with hadronic systems under extreme conditions where 
the interactions are mediated by the strong force (see 
Ref. (46J), but no experimental signature(s) of its real- 
ization has been proposed so far. Atomic systems of- 
fer a novel setting for studying the LOFF phase under 
conditions that are more favorable than those in solids 
(absence of the lattice effects, access to the momentum 
distribution in the system through time-of-flight experi- 
ments) mi H3 HE HI US- Atomic systems also offer the 
possibility of novel realizations of the LOFF phase which, 
for example, invoke P-wave anisotropic interactions |4l| 
or finite-size systems [4J, |45j ■ 

While the LOFF and DFS phases break global space 



symmetries, the way they do so is fundamentally dif- 
ferent: the LOFF phase breaks both the rotational and 
translational symmetries due to the finite momentum of 
the condensate and irrespective of the form of the lat- 
tice structure; the DFS phase breaks only the rotational 
symmetry from 0(3) down to 0(2). We shall compare be- 
low the realizations of the DFS and LOFF phases in the 
6 Li gas assuming that the order parameter in the LOFF 
phase has a simple plane-wave form; in addition, we shall 
describe an experimental signature of the DFS phase that 
can be established in time-of-flight experiments and that 
would allow one to distinguish the DFS phase from the 
competing phases. A brief account of this argument was 
given earlier in Ref. 47] . Yet another possible alternative 
is the phase separation of the superconducting and nor- 
mal phases in real space, such that the superconducting 
phase contains particles with the same chemical poten- 
tials - i.e. is symmetric - while the normal phase remains 
asymmetric; see Refs. .48, ,4Sj|£3- A comparison of the 
heterogeneous phase with the alternatives would require 
knowledge of the poorly known surface tension between 
superconducting and normal phases, and will not be at- 
tempted here. 

This paper is organized as follows. Starting from 
Dyson's equations, we derive in Sect. II the dispersion 
relations of elementary excitations in an asymmetrical su- 
pcrfluid. Sect. Ill discusses the gap equation and its reg- 
ularization for contact interactions. Sect. IV is devoted 
to modifications of the single particle spectra needed to 
describe the LOFF and DFS phases. Numerical results 
for the case of 6 Li gas are presented in Sect. V. Sect. VI 
discusses the possibility of observing the DFS phase in 
time-of-flight experiments. 



II. FORMALISM 



Consider a uniform gas of fermionic atoms with two hyperfine states, which we assign labels f and J. (these states 
equivalently can be thought of as pseudospins.). The model Hamiltonian that describes our system is 

^ = ^ E / d 3 x\/^ a (x)V^ a (x) -Y,f d3x [ d 3 x'i>Ux)i>l(x')V(x, x')fo(x')$ a (x), (2) 

a a/3 

where ^ a {x) and ip a {x) are the creation and annihilation operators of a state at the space point specified by the 
position vector x and pseudospin a{~], J.), m is the atom's bare mass, and V(x,x') is the two-body potential (here 
and below the volume f2 = 1). The one-body propagator in the superfluid state is a 2 x 2 matrix in Gor'kov space, 

r( , ( G aP {x,x') F aP {x,x')\ _( -i(T^ a {x)^{x')) (V»„(x)V/j(^)> \ r , 
~ [XlX> \-Fl fi {x,x>) G a0 {x,x') ) \ (^(x)4(x')) -i(f M x)4(x')) ) ' {6) 

where x — (x,t), the Greek indices stand for the pseudospin, T and T are the time-ordering and time-antiordering 
symbols, and G a p(x,x') and F^(x,x') are the normal and anomalous propagators The propagators above are 
assumed to be ordered on the Schwinger-Keldysh real-time contour 52], which permits a finite-temperature treatment 
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of the problem. To keep our presentation concise, we shall not write down the explicit form of the correlation functions 
on the contour; our final expressions are written for retarded correlation functions. The 2x2 matrix Green's function 
satisfies the familiar Dyson equation 



G_ a p{x, x') — G®,p(x, x') + ^2 Jd 4 x"d 4 x'"G^ ll (x 1 x"')^ lS (x"' , x")G g ^(x" ,x'), 

7,5 



(4) 



where the free propagators G^p(x, x') are diagonal in Gor'kov space; the underline indicates that the propagators 
and self-energies are matrices in this space. The components of the Fourier transform of Eq. (0J with respect to the 
difference of the space arguments of the two-point correlation functions are 



G a p{p) = G a/3(p) + G 0ai (p) T, lS (p)G S p{p) + A lS {p)F} (p) 



(5) 
(0) 



where p is the four- momentum, Go a p (p) is the free nor- 
mal propagator, and and A a p(p) are the normal 
and anomalous self-energies; summation over repeated 
indices is understood. The Dyson equations for the com- 
ponents G a p(p) and F a p{p) follow from Eqs. (0 and © 
through the time-reversal operation. 

Below, we shall assume that the interactions conserve 
spin- i.e., G a p(p) = 5 a pG{p) and T, a/3 (p) = <5 Q/3 £(p) - 
and concentrate on the pairing in the state of zero total 
spin and orbital angular momentum S = L = 0. Thus 
the anomalous propagators and self-energies must be an- 
tisymmetric with respect to the spin indices, 



fL(p) = g a0 F\ P ), A* (p) = g aP A\ P ), 



(7) 



where g a p = ia y is the spin matrix, with a y being the 
second component of the vector of Pauli matrices. It 
is convenient to rewrite Eqs. JSJ and ® in terms of 
auxiliary Green's functions, which describe the unpaired 
state of the system: 

G&(p) = G 0af) (p) + G^(p)^ s (p)G^ f3 (p). (8) 

The formal solution of Eq. (JSJ for the retarded propaga- 
tors in terms of the self-energy E(p) is 



'±0J - £ pcr - Scr(±p)] 



d a p[±(u> + irj) 



0(ImS CT (±p)), 



(9) 



where £ p(7 = p 2 /2m a — [i a is the energy of a free spin- 
<7 particle relative to the chemical potential fj, a ; the sec- 
ond line follows in the quasiparticle approximation, which 
keeps the leading-order term in the expansion of the prop- 
agator with respect to the small imaginary part of the 
self-energy. The quasiparticle dispersion relation in the 
normal state is then given by e pa — £ p(T + ReS CT (±p). 
Combining Eqs. JSJ and © with Eq. JHJl one finds 

G aP {p) = G^( P ) [5 1(3 + A lS (p)F^(p)] , (10) 
fUp) = G^{-p)A\ s {p)G 50 { P )- (11) 



these equations are easily solved to obtain the propaga- 
tors 



G U (P) = 



u> + E s ±E A 
{lo - E A ) 2 -E 2 -A 2 ' 



fHp) = - 



At 



(w - E A ) 2 - E% - A 2 • 



(12) 
(13) 



where E s = (e pT + e p j.)/2 and E A = (e pT - e p \)/2 are, 
respectively, the parts of the spectrum which are sym- 
metric and antisymmetric under time-reversal operation. 
The poles of these propagators (which must be identical) 
define the dispersion relation of the quasiparticles in the 
paired state: 



^1/2 = E A ± 



El 



(14) 



Note that the quasiparticle spectrum is twofold split due 
to the asymmetry in the number of the spin-up and spin- 
down atoms. In the symmetric limit (E A — 0) we recover 
the ordinary BCS dispersion relation. 



III. GAP EQUATION 

To obtain a closed set of equations we need to specify 
the approximation to the anomalous self-energy. In the 
mean-field (BCS) approximation, 



At(p) = i / V{p,p')F\p') 



d 4 p' 
(2tt)4 



i / V(p,p')G? (-p')AV)G T (p') 



d 4 p' 

(15) 

where V(p, p') is the Fourier transform of the two-body 
interaction V(x,x') which is responsible for the pairing. 
The pairing interaction can be renormalized in such a way 
that integration in Eq. I|15|) is restricted to the vicinity 
of the Fermi surface. This permits us to approximate 
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smooth functions of momentum by their value at the 
Fermi momentum. We introduce a momentum renor- 
malization scale A such that A <c e\ <C min[e^|, £fj.], 



where £f\i are the Fermi energies of spin- up and down 
species. Then, 



At(p) = i J U{p,p')G^{-p')^{p')G^ P ')e(K-\p'\)^- A , (16) 

U(p,p') = V(p,p') + i J V (p,p")Gf(~p")Gf(p")U(p",p')9(\p"\-A)0^. (17) 

In the second equation the full propagator G|(p) was replaced by its counterpart in the unpaired state, G^ (p). By 
construction, the on-shell integration is carried over momenta much larger than the Fermi momentum for which the 
quasiparticle spectrum is unaffected by the pairing gap. Thus Eq. 1(17(1 decouples from Eq. (|16fl . while the integration 
in Eq. 1(16(1 is now constrained to the vicinity of the Fermi surface. This permits us to simplify the problem by 
expanding the pairing interaction in spherical harmonics with respect to the angle formed by the momenta lying on 
the Fermi surface(s) and, thus, to reduce the pairing interaction to a function of a single angle. While this is useful 
for finite-range interactions, in the case of zero-range interactions the potential V(p,p') needs to be eliminated from 
Eq. Ijl 7(1 in favor of the scattering T matrix, which obeys the integral equation 

T(p, P ') = V{p 1 p')+ l J^ V (p,p")Gf(-p")Gf(p"W,p')j^. (18) 

Combining Eqs. 1 (17(1 and ((18(1 one finds a regular integral equation defining the pairing force 

U(p,p') = T{p,p')- l j\{p,p'')G N i ^p' l )G^{p' l )T{p'\p')^. (19) 



In the dilute limit of interest, partial waves higher than 
the s-wave can be neglected, and the interaction is solely 
determined by the s wave scattering length as < 0, as 
T(p,p') = T = 47r|a s |/m. The solution of Eq. 
with this interaction is straightforward, and one obtains 
iorU(p,p') = U 



The densities of the species are given by 

d 3 p'dui' 



P\li = \ Im G TT / u (p',o;)/(w)- 



(2tt)< 



{U(p') 2 [f(L0 Vl ) - f(-LO m )] + /(-W; /T )> 



p' 2 dp' 
(22) 



U = T 



l-/Gf(-p)Gf(p)fl(A-N)^4 



.(20) 



For the zero-range interaction above the gap equation 
takes the form 



At(p) = U J ImFtfcVJjVMA-lp'l) 



A 



Uo 



2 Jo y/EfTE? 



</(wi) - /(wa)> 



,^ d 3 p'duj' 

P ,2 dp' 



(27T) 



(21) 



where (...) stands for the angle average and the second 
line follows in the quasiparticle approximation [i. e., by 
retaining only the pole part of the propagator F'(p)]. 



where u{p) 2 = 1/2 + E S /(2 V / E 2 + A 2 ) is the familiar 
Bogolyubov amplitude and the second line of Eq. IL'21) 
is obtained in the quasiparticle approximation. The nor- 
mal self-energy in the T-matrix (ladder) approximation 
is defined as 



T aa ,{p,p'-e p + u')G a ,{p') 



(2tt) 4 



(23) 



where the dependence of the self-energy on the center- 
of-mass momentum is suppressed. If the T matrix is 
approximated as above by a real constant To which ac- 
counts for interactions between different species, the nor- 
mal self-energy is momentum independent, purely real, 
and is given by 



s T/i - T oPl/^■ 



(24) 
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This constant shift in quasiparticle energy can be ab- 
sorbed in the chemical potential by defining A**/j, = 
fiyi — Sjrt. It is straightforward to extended this re- 
sult to the case where the T matrix contains components 
which act among the same species or includes contribu- 
tions from partial waves with L > 0. Equations Q2()jl. 
(|21|1 , 11' I'D , and (|24[1 are the fundamental equations of the 
mean-field theory for systems with unequal spin pairing 
that are interacting via zero-range forces. They include 
implicitly the effects of a finite momentum of Cooper 
pairs and the topology of the Fermi surfaces to which we 
shall turn in the following section. 



IV. FINITE PAIR MOMENTUM AND 
DEFORMED FERMI SURFACES 

While the BCS ground state assumes that fermions 
bound in a Cooper pair have equal and opposite momenta 
(and spins), for fermionic systems with and unequal num- 
ber of spin-up and -down particles this is not always true. 
Larkin and Ovchinnikov |39j and independently Fulde 
and Ferrell 01 observed that pairing is possible among 
pairs which have a finite total momentum with respect 
to some fixed reference frame. The finite momentum P 
changes the quasiparticle spectrum of the paired state. 
To see this note that the auxiliary propagator writ- 
ten in the center-of-mass frame reads 



1 ( P + 
u> — ± p 

2m V 2 



I'Vl + iV 



,(25) 



and, therefore, the symmetric and antisymmetric (under 
time-reversal) parts of the quasiparticle spectrum are 



E S = - 



1 P 1 



4m 



w 1 f Pp *^ * 



/' - - /' j (26) 
(27) 



The results of Sec. IIIII remain valid with the above re- 
definitions of Es and Ea- Note that the quantities of in- 
terest, in particular the gap, now depend parametrically 
on the total momentum. As a consequence, the twofold 
splitting of the spectrum i|14|) does not vanish in the limit 
of an equal number of spin-up and -down particles (i.e., 
when /j,* = p,*). While such a state lowers the energy of 
the system with respect to the unpaired state, it is still 
unstable with respect to the ordinary BCS ground state. 

We now turn to the deformations of the Fermi sur- 
faces. The two Fermi surfaces for spin-up and -down 
particles arc defined in momentum space by the equa- 



tions £ 



T/J 



£p/2± P — Vyi 



= 0. When the states 



are filled isotropically within a sphere, the chemical po- 
tentials are related to the Fermi momentum pf,g as 
p,* = pp a ./2m (for the sake of argument we assume here 
that the temperature is zero). To describe the defor- 
mations of Fermi surfaces from their spherical shape we 



expand the quasiparticle spectrum in spherical harmon- 
ics e a — J2i £ io-Eii x )i where x is the cosine of the angle 
formed by the particle momentum and a randomly cho- 
sen symmetry-breaking axis and Pi(x) are the Legendre 
polynomials. The I — 1 terms break the translational 
symmetry by shifting the Fermi surfaces without deform- 
ing them; these terms are ignored below. Truncating the 
expansion at second order (I = 2), we rewrite the spec- 
trum in a form equivalent to the above one |38| 



-T/l 



ep± p -M*/i (i + ^T/l^ 2 ) 



(28) 



where the parameters r\ a describe the quadrupole de- 
formation of the Fermi surfaces. It is convenient to 
work with the symmetrized 5 = [rj^ + T)i)/2 and anti- 
symmetrized Se = (r)-\ — f)i)/2 combinations of t^j.. Be- 
low we shall assume S = and consider two limiting cases 
Se ^ and P — (the phase referred as the deformed 
Fermi surface superfluidity (DFS) phase) and Se = and 
P^O (the plane- wave LOFF phase). 



V. RESULTS 

Consider a trap loaded with 6 Li atoms and assume 
that the net number of atoms in the trap is fixed while 
the system is maintained at constant temperature. As- 
sume further that the number of atoms corresponds to a 
Fermi temperature Tp — ep/ks = 942 nK, which in the 
uniform, symmetric case at T = would translate into a 
Fermi momentum of the system kp ss 4.83 x 10 4 cm -1 
and a density p — 3.8 x 10 12 cm~ 3 . We shall work below 
at constant temperature T = 10 nK (<C Tp); i.e., the 
system is in a hig hly degenerate regime. (In the experi- 
ments of Ref. |53j , such a Fermi temperature corresponds 
to about 4xl0 5 atoms in a single hyperfme component.) 
Experiments control the partial densities of atoms in two 
different hyperfme states |f ) = \F = 3/2, nip = 3/2} and 
||) = |3/2, 1/2), e. g. by transferring atoms from one 
state to the other using ~76 MHz rf pulses [54[. Since 
the free-space triplet scattering length for 6 Li atoms in 
these hyperfme states is as — — 2160as, the system is in 
the weakly coupled regime fcjr |ag | <^7r/2. 

The pairing gaps of the LOFF and DFS phases com- 
puted from the coupled equations (|2*T)l and are 
shown as a function of asymmetry parameter a = (pf — 
Pi)/(PT+Pl) f° r different values of the total momentum P 
and deformation 8e in Fig. ^ Without loss of generality 
the density asymmetry is constrained to positive values; 
i.e., we assume p-^ > py The positive values of 5e corre- 
spond to a prolate (cigar like) deformation of the major- 
ity and oblate (pancake like) deformation of the minority 
population's Fermi spheres; we shall confine ourselves to 
the case where Se > since we have checked that it is 
the one corresponding to the largest value of the pairing 
energy. 

To elucidate the dominance of the phases with broken- 
space symmetries over the asymmetric BCS state con- 
sider the modifications in the single-particle spectra im- 
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FIG. 1: (Color online) The dependence of the pairing gaps 
in the LOFF phase (upper panel) and the DFS phase (lower 
panel) on the asymmetry parameter for several values of the 
the total momentum P/Icf and deformation parameter 8e 
which are indicated in the panels. 
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FIG. 2: The dependence of the quasiparticle spectra of two 
hyperfine states oji and u>2 on the momentum for a = = Se 
(solid line), a — 0.05 and Se = (long-dashed lines), a = 0.05, 
and Se = 0.1, x = (short-dashed lines), and x = 1 (dash- 
dotted lines). The Fermi momentum Rf = 4.83 x 10 4 cm -1 
is indicated by the vertical line. 



plied by these phases. In the asymmetric BCS state, 
the antisymmetric part of the quasiparticle spectrum 
(|14fl . Ea, acts in the gap equation (|21|l to reduce the 
phase-space coherence between the quasiparticles that 
pair (when Ea — the BCS limit is recovered with equal 
occupations for both particles and perfectly matching 
Fermi surfaces). This blocking effect is responsible for 
the reduction of the gap with increasing asymmetry and 
its disappearance above a ~ 0.055. 

When the pairs move with a finite total momentum or 
the Fermi surfaces are deformed, the antisymmetric part 
of the spectrum, Ea, is modulated with the cosine of the 
polar angle x (in the frame where the z axis is along the 
symmetry-breaking axis). In the case of the plane- wave 
LOFF phase Ea oc x while in the DFS phase Ea oc x 2 . 
This variation acts to restore the phase-space coherence 
for some values of x at the cost of even lesser (than in 
the BCS phase) coherence for the remaining values. The 
effect can be seen explicitely in Fig. which compares 
the quasiparticle spectra uj\ and 0J2 in the DFS phase 
for combinations of a and Se for two orthogonal direc- 
tions with the BCS result for the symmetric case (black) . 
Along the symmetry-breaking axis the energy separation 
between the quasiparticle spectra is considerably smaller 
than in the asymmetric BCS state; in the orthogonal di- 
rection, the opposite is the case. Compared to the asym- 
metric BCS state the phase-space overlap between pairs 
is increased in the first case and decreased in the second. 
The net result, displayed in Fig- El is the increase in the 
value of the critical asymmetry a c at which the super- 
fluidity vanishes. At large asymmetries the DFS phase 



exhibits the reentrance effect: the pairing exists only for 
the deformed state between the lower (a cr i) and upper 
(«cr2) critical deformations; see Fig. |3| This figure also 
shows the largest possible asymmetry supported by the 
LOFF for a given total momentum P; this asymmetry is 
maximal at P/kp ~ 0.075. 

An important feature of the spectrum of the asym- 
metric BCS state (a _/0, Se = = P) is its gapless 
nature [5^|5^|52i|5l,|5jj- i.e., the existence of nodes for 
one (or both) branches of the spectra (cf. the gapped 
BCS spectrum also shown in Fig.|2J. Gapless excitations 
affect the dynamical properties of the superfluid state 
such as the transport and collective modes and lead to 
a number of peculiarities in the thermodynamics of this 
state. This feature clearly remains intact for the phases 
with broken-space symmetries. As seen in Fig. [21 the 
spectrum of the DFS phase covers a range bounded by 
the curves with x — and x — 1 and features nodes at 
which the quasiparticles can be excited by an infinitely 
small external perturbation. The macrophysical manifes- 
tations of the LOFF and DFS phases such as the response 
to density perturbations or electromagnetic probes and 
the thermodynamic functions (heat capacity, etc.) would 
differ from the ordinary BCS phase due to the nodes and 
anisotropy of their spectrum. Such an anisotropy can 
be used to discriminate phases with broken-space sym- 
metries in the time-of- flight experiments (see Sect. IVl|l. 
Moreover, the phases with broken-space symmetries fea- 
ture a larger number of Goldstone modes than the asym- 
metric BCS phase because of the breaking of additional 
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FIG. 3: The dependence of the critical asymmetry a c of 
the transition from the superfiuid to the normal state on the 
total momentum in the LOFF phase and the deformation 
parameter in the DFS phase. The model parameters are as 
in Fig. 1. 



FIG. 4: (Color online) The dependence of the free energy 
of the LOFF (upper panel) and the plane-wave DFS phase 
(lower panel) on the asymmetry parameter for several values 
of the deformation parameter Se and the total momentum 
P/fc_F which are indicated in the panels. 



global space symmetries. 

Which phase is the true ground state at a give den- 
sity asymmetry is decided by a comparison of the free 
energy difference AF = F$ — -F/v of these phases which 
are displayed in Fig. 0] The LOFF phase is preferred 
to the normal and homogeneous BCS phases in a nar- 
row window of asymmetries, 0.04 < a < 0.057, and for 
the total momentum of the pairs, P/kp ~ 0.05. This is 
consistent with the results obtained in the scheme where 
the density asymmetry is described in terms of the dif- 
ference in the chemical potentials of the species <5/i [the 
critical value for the BCS phase is S^f cs = 0.707A(0), 
while for the LOFF phase S[i L c OFF = 0.755A(0), where 
A(0) ee A(Sfi = 0)] ^1113. Note that while there is a 
nontrivial solution to the gap equation for P/Uf > 0.01, 
the gain in the pairing energy is less than the loss in ki- 
netic energy due to the motion of the condensate and the 
net energy of the LOFF phase is larger than that of the 
asymmetric BCS phase. However, the pairing energy of 
the LOFF phase can be increased by choosing a more 
complex form of the order parameter - e. g., by keeping 
a larger number of terms in the expansion of the order 
parameter in the Fourier series. 

The DFS phase is the ground state of the system 
(i.e., it has a lower free energy than the normal, BCS, 
and LOFF phases) in a wider range of asymmetries, 
0.03 < a < 0.06, for the deformation parameters in the 
range 0.12 < Se < 0.16. For even larger deformations the 
gain in the pairing energy does not compensate the loss 
in the kinetic energy due to the stretching of the Fermi 



surface into ellipsoidal form. Note that the free energy 
is also affected by the reentrance effect (i. e. restoration 
of pairing correlations as the asymmetry is increased). 
However, while beyond the reentrance point the DFS 
phase becomes preferable to the unpaired state, its free 
energy is still larger than in the homogenous asymmetric 
BCS state. Only at the large asymmetries quoted above 
does it become the ground state of the system. 

To summarize, the coherence is restored and the 
strength of pair correlations is increased in the LOFF 
phase due to the finite momentum of the Cooper pairs. 
In the DFS phase the same is achieve by stretching the 
spherical Fermi surfaces into ellipsoids. The fundamen- 
tal difference between these phases is that the transla- 
tional symmetry remains intact for the DFS phase, which 
breaks only the rotational symmetry, while the LOFF 
phase breaks both symmetries. Quantitatively, the maxi- 
mal value of the gap and the absolute value of the ground- 
state free energy is larger in the DFS phase than in the 
LOFF phase for asymmetries larger than a ~ 0.04. For 
these asymmetries both phases are favorable over the ho- 
mogeneous BCS phase. However, one should keep in 
mind that the LOFF phase admits a variety of lattice 
forms and the plane-wave structure need not be the most 
favored one |46J. 
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FIG. 5: (Color online) The dependence of the occupation 
probabilities of two hyperfine states on the momentum. The 
Fermi momentum fci? = 4.83 x 10 4 cm -1 is indicated by the 
vertical line. The labeling of the lines is as follows: a = 
= 8e (solid line), a = 0.05 and 8e = (long-dashed lines), 
a — 0.05, 5e — 0.1, x — (short-dashed lines), and x = 1 
(dashed-dotted lines). The bell-shaped light curves show the 
anisotropy - the difference between the x = 1 and x — 
occupation numbers - for a = 0.05, 8e = 0.1. The remaining 
parameters are as in Fig. 1. 



VI. DETECTING THE DFS PHASE IN 
EXPERIMENTS 

The large number of side effects in solid systems (such 
as the presence of the ionic lattice or defects) have pre- 
cluded up to now a clear detection of a superfluid phase 
with broken-space symmetries in these systems. The sit- 
uation can be much more favorable in ultracold gases 
due to the high control over the experimental conditions. 
Recently, Mizushima et al. proposed a way to detect the 
LOFF phase in these systems by i mag ing the density pro- 
files by local spin measurements |44j . Below we propose 
simple method to detect the DFS phase in a Fermi system 
with rotational symmetry - e.g., a system which is cither 
homogeneous or in a spherical trap. Indeed, experimen- 
tal evidence for the phases with broken-space symmetries 
can be obtained from studies of their momentum distri- 
butions which, unlike in the homogeneous phase, must 
be anisotropic in space. Figure [5] shows the occupation 
numbers in the BCS, the asymmetric BCS, and the DFS 
phases; varying the cosine of the polar angle x covers a 
range of occupation probabilities which includes the un- 
deformed asymmetric state. The bell-shaped curves show 



the angular polarization of the occupation numbers in the 
DFS phase defined as 8n a = |n CT (a; = 1) — n a (x = 0)|. 
The maximal anisotropy in the occupation probabilities 
of particles along and orthogonal to the symmetry break- 
ing axis is about 90%. Thus, a direct way to detect the 
DFS phase is the measurement of the anisotropy in the 
momentum distribution of the trapped atoms. Such a 
measurement can be realized by the time-of-flight tech- 
nique [IHSSHIiSlil. This method uses the fact 

that after releasing the trap, the atoms fly out freely and 
an image of their spatial distribution taken after some 
time of flight provides information on their momentum 
distribution when confined inside the trap. Assuming 
that the system was in the deformed superfluid state 
one would detect a mean momentum of particles of type 
1 (majority) in the direction of symmetry breaking by 
about 90% larger than that of particles of type 2 (mi- 
nority) in the same direction. Therefore, the presence of 
this anisotropy in the detected momentum distributions 
would be evidence for a deformed superfluid state be- 
ing the ground-state of the system, as deformation alone 
(i, e. without pairing) would not lower the energy so as 
to produce a deformed nonsuperfluid ground-state. Note 
that this argument is equally valid for a homogeneous 
system or for an atomic gas in a spherical trap, where no 
preferred direction is introduced by the trapping poten- 
tial. For a nonsphcrical trap, the momentum distribu- 
tions of both species are expected to be deformed in the 
same way; therefore, the detection of an anisotropy in 
momentum distributions superimposed on the deforma- 
tion induced by the deformed trap could be a signature 
of the DFS phase. Case studies for certain forms of de- 
formed trapped (e.g., in the local density approximation) 
would be necessary to quantify the effect of combined de- 
formation due to the trap and pairing. 

The direction of spontaneous symmetry-breaking (in k 
space and, therefore, also in real space) is chosen by the 
system randomly and needs to be located in an experi- 
ment to obtain maximum anisotropy. A clear distinction 
between the DFS and LOFF phases can be achieved in 
time-of-flight experiments, since the latter predicts pe- 
riodic momentum distributions, independent of the de- 
tailed spatial structure of the phase that nucleates in the 
ground-state. 
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